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Abstract 

The method for the recursive calculation of the effective potential is 
applied successfully in case of weak coupling limit (g tend to zero) to a 
multidimensional complex cubic potential. In strong-coupling limit (g tend 
to infinity), the result is resumed using the variational perturbation theory 
(VPT). It is found that the convergence of VPT-results approaches those 
expected. 
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1 Introduction 

Quantum dynamics allows to avoid the operator formalism for the ben- 
efit of infinite products of integrals, called Path-Integrals [1,2]. Unfortu- 
nately, most Path-Integrals can not be performed exactly; therefore, many 
different approximation procedures are developed in order to deal with 
non-analytically solvable systems. 

The most commonly used one is known as perturbation theory [3-9] . It 
is based upon the expansion of some physical quantity, e.g. the ground- 
state energy of particle in some potential, into a power series of coupling 
constant, the results obtained in weak-coupling limit seem to converge 
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to the exact result for low orders, so that the divergence of that series 
becomes important when the expansion is driven to the higher-order [3- 
12]. Therefore, it is necessary to find means for treating the divergent 
perturbation series, including the strong-coupling limit. Thus, we apply 
the resummation scheme often based on introducing an artificial parameter 
trick, namely Kleinert's square-root [2,13-15]. 

The variational perturbation theory (VPT) [2,7-9] has been developed 
thanks to the variational approach due to R. P. Feynman and H. Kleinert 
[6], and has been extended to an efficient non-perturbation approxima- 
tion. VPT allows the conversion of the divergent weak-coupling into the 
convergent strong-coupling expansion. 

According to some studies, a more physical alternative axiom called 
space-time reflection symmetry (VT— symmetry) allows for possibility of 
cubic complex non-Hermitian Hamiltonian but still leads to a consistant 
theory of quantum mechanics. It is shown that if the VT — symmetry of 
a Hamiltonian H is not broken, then the accompanying spectrum remains 
real and positive. 

The space reflection operator V is a linear operator with the property 
V 2 = 1 and has the effects p — > —p, and x — > —x, while the time reflection 
T is an antilinear operator with the property T 2 = 1 and has the effects 
p — > —p, x — > x, and i — > —%. 

Recently, the convergence of VPT has been tested successfully for the 
ground-state energy of the one-dimensional complex cubic coupled oscilla- 
tor [11,13,15] 

V (x) = ^-co 2 x 2 + igx 3 . (1) 

The purpose of this paper is to examine how, in a first approach, the VPT 
can be applied to the generalized 2D and 3D-complex cubic potentials 

M 

V(x,y) = —u 2 (x 2 + y 2 )+igxy 2 . (2) 

V (x, y, z) = yw 2 (x 2 + y 2 + z 2 ) + igxyz. (3) 

by resuming the weak-coupling series [10,12] of the ground-state energy via 
VPT. 



2 



In the strong-coupling limit, the potentials (2) and (3) are reduced 
to potentials without a harmonic term. It turns out that the rate of con- 
vergence is not satisfactory. Therefore, combining the effective potential 
[13-15] to the VPT permits the improvement of the rate of convergence. 
This later should be performed to higher orders, and since mathematical 
computation become cumbersome, we have stopped our calculation at the 
second order in h. 



2 Perturbation theory : Feynman diagrams 



In this section, we derive the weak-coupling coefficients for the ground-state 
energy of the potentials (2) and (3) using the Feynman diagrammatical 
expansion. The partition functions 




will be calculated perturbatively by an expansion in the coupling constant 
g, here we have used x a (resp. x a ) instead of x a (r) (resp. x a (r) ), and 
dot {< ' * refers to derivatives of x a (r) with respect to r. These expressions 
can be converted, using the Taylor expansion of the exponential function, 
into the form 
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Introducing in both cases the notation, namely expectation values, 
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lead, up to the second order in g, to the following expressions of the par- 
tition functions 
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where Z^ are the partition functions of the harmonic oscillator 



Z- = Z- = Z«=i S inh-'^. 
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Expectation values in (9) and (10) can be performed by applying general- 
ized Wick's rules 



• each expectation values of a product of 2-paths defines the propagator 

(Xi (ti) Xj (t 2 )) u ^ 0; for i = j. 
(Xi (ri) Xj (t 2 )) u = 0; for i ^ j. 



G w (ti,t 2 ) = 



(12) 



• each expectation values of a product of n-paths reads, taking into 
account the first rule (12) 

(xi (n) Xj (r 2 ) • • • x q (r„)) w = G w (n, r 2 ) (x fc (r 3 ) x z (r 4 ) • • • a, (r n )) w 

+ G w (n, r 3 ) (xj (r 2 ) • • • (r n )) w H 

+ C w (ri, r„) (xj (r 2 ) • • • rr p (r„_i)) w , 

(13) 

By applying the Wick's rules and knowing that the free energy reads 
as F = —ksT \nZ, we derive the perturbation series at low temperatures 
(T — > 0) for the ground-state energy of the potentials (2) and (3) up to the 
fourth-order by the connected Feynman diagrams 



E- 
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(14) 



0{g 6 ). (15) 



For evaluating the connected Feynman diagrams, the following Feynman 
rules are introduced [2,14]: 



5 



• for propagator 



h 

(T a ,r b ) = 7^ exp [-u \r a - r b \] . (16) 



for vertices 

-6ig M 

~* h 



dr a . (17) 



Applying the Feynman Rules to (4) and (5), we obtain the analytical 
expressions for the ground-state energy [10,12] 



F (o) _ ^ , 5 9 2 h 2 223 gW 6 



E. 



3D = ^ + 
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+ O (g 6 ) . 



(19) 



3 Bender- Wu perturbation theory : Recur- 
sion relations 

We derive here recursion relations of the perturbation coefficients of the 
ground-state energy for the potentials (2) and (3). These recursion rela- 
tions are obtained from the corresponding Schrodinger equations 



2M 



fay) + 



h 2 



2M 



2 

Muj 2 



(x 2 + y 2 ) + igxy 2 



(x 2 + y 2 + z 2 ) + igxyz 



if) fay) = E 2D i) (x,y). 

(20) 

il>(x,y,z) = E ZD i\) fay,z) 



(21) 

In the perturbation theory, the ground-state wave functions are ex- 
panded in the form [3-5,12] 



ifj (x, y) = j\f exp 



Muj 



(x 2 + y 2 ) + <f>(x,y) 



(22) 
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ip(x,y,z) = A/" exp 



~2/T 



(x 2 + y 2 + z 2 ) +<j>(x,y,z) 



(23) 



The wave functions <f> (x a ) , where a = 1, 2, 3, will be expanded in pow- 
ers of the coupling constant g 



(24) 



k=i 



The ground-state energy corresponding to (2) and (3) can be expanded 
in powers of the coupling constant g 



E 2D = huj + Y,3 k 4 



k(2D) 
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kj3D) 
k 



(25) 



(26) 
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Inserting (24), (25) and (26) into (22) and (23), we obtain a differential 
equation taking into account natural units {h = M = uj = 1) 

e[ 2D) = ~ V^A ( x a) + x d x (j) k (x a ) + y d y <f> k (x a ) + igxy 2 5 kA 

j fc-1 

~2 [ d ^k-i ( x a) d x </>i (x a ) + d y (p k ^ (x a ) dyfa (x a )] . (27) 
i=i 
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where we have used the abbreviation ct = 



_d_ 

8x r 



The 6 h ( expanded in powers of the coordinates as [3-5,12] 
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^(^) = EE4!^v m (29) 

<f) k (X, l/,*)=EEE otl^y^- (30) 
j=0 m=0 n=0 

where and are non-symmetrical coefficients and can be real 

anc?/ or imaginary. 

By inserting (29) and (30) into (27) and (28), we obtain in second order 
the following coefficients 
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For k — 1, 
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(2) _ (2) _ (2) _ (2) _ (2) _ (2) 
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= ~, all others af^ n = 0. 



72' 

(34) 

For k > 3, we find the recursion relations and the energy correction 
coefficients 
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Table 1 and Table 2 show the coefficients up to the 10th order [10,12]; 
one sees that the results obtained by the Bender- Wu method are in agree- 
ment with those obtained by the connected fourth order Feynman diagrams 
(see (18) and (19)). 

Table 1 

Weak-coupling coefficients for the 2-dimensional potential up to the 10th 
order, k represents the order of expansion and coefficients of the energy 
corrections. 
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2 


3 


4 
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ffc 





5 

24 





223 
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155520 





346266143 
111974100 





2360833242959 
141087744000 
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Table 2 

Weak-coupling coefficients for the 3-dimensional potential up to the 10th 
order, k represents the order of expansion and coefficients of the energy 
corrections. 



k 


1 


2 


3 


4 


5 







i 

24 





5 

576 
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7 
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10 
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5069 





2441189 





8034211571 


622080 


289598400 


752467968000 



4 Resummation procedure : Strong-coupling 
limit 

In this section, we are interested to resume the perturbation series (18) and 
(19) to the strong-coupling limit. To this end, substituting the coordinate 
x a , where a = 1, 2, 3, by 



x a > g ^ x a , 



the Schrodinger equations (20) and (21) become 



h 2 



2M 



(x,y)+ 



M u< 



9 ' 5 (x 2 + y 2 ) + ixy 2 



(39) 

ifj(x,y) = g-^ 5 E 2D ij (x,y) 
(40) 



h 2 



2M 



M u' 



g~ 4/5 (x 2 + y 2 + z 2 ) + ixyz 



tfj(x,y,z) 



g-^E 3D i;(x,y, 



(41) 



The wave functions and the energy will be expanded in powers of the 
coupling constant g and yield [13,15] 

i) (x a ) = % (x a ) + <T 4/5 ^i (x a ) + <T 8/5 ^2 M + ■■■ ■ (42) 
E D = g 2 / 5 (e + g-V 5 e 1 + g-*/ 5 e 2 + ■■■). (43) 
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The strong-coupling coefficients of the ground-state for the potentials 
(2) and (3) can be obtained by resuming the weak-coupling series [10,12] 
obtained in the last two sections. Up to the fourth order, the weak-coupling 
series read 

„foi 5 g 2 h 2 223 g 4 h 3 „ , fi . 

(0) 3 1 g 2 h 2 7 g*H 3 6 

4i = - 576 W + ° k ) • (45) 

The power behavior of taken in (43) is independent of the order 
iV and solely the coefficients depend on N 

E (m = g y 5 + g -4,s € (N) + g -s /5e im + . . . ) (46) 

In order to perform a resummation, an artificial parameter trick is 
introduced often called Kleinert's square-root [2,13-15] 

u -> fi^T+^V, (47) 

with 

(48) 

Therefore, inserting (47) into (44) and (45) and re-expanding in cou- 
pling constant g 2 to the first order taking into account (48), we obtain 



sS( S ,^? + ^ + ^ + 0( S <). (49) 



„(i") / 3KI 3hio 2 oh 2 „ / a\ / 

fo, fl) = — + — + g 2 — + O ( 5 4 ) . (50) 

Performing the first derivative with respect to the variational parameter 
Q and considering its strong-coupling behavior [13,15] 

Q = g 2 ' 5 (fi + (T^i + c7- 8/5 fi 2 + • • • ) (51) 

we get 
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• 2-D 



<5h „ lo 2 J 3 u A J 27 



il o = \— ; Hi = — (/— ; tt 2 = —{ — . (52) 

"3 5 V 5ft ' 25 V 125ft 3 V ; 



• 3-D 

^2ft u; 2 J 9 „ u; 4 J 3 



° V 9 ' 5 V 2ft ' 25 V 8ft 3 v ; 

Inserting these results into (49) and (50) yield the strong-coupling co- 
efficients of the ground-state energy ejjf^ 

• 2-D 



a. = 0.69222/^; «f» = ^{/^; 4" " * /3S2 



8 V 3 1 2 V 5 z 20 V 5 

(54) 



• 3-D 



(55) 

However, it turns out that the numerical values of the leading strong- 
coupling coefficients obtained in (54) and (55) are lower than those calcu- 
lated in weak-coupling limit (see (18) and (19)). Thus, the convergence of 
the VPT-results is less satisfactory. Following [13], we must rely upon that 
the farfetched values must be higher than those calculated in weak-coupling 
limit. It is what we are going to show in section 5. 



5 D-dimensional complex effective potential 

For any given potential, the corresponding effective potential is obtained 
by performing a Legendre transformation of the ground-state energy in 
the special case of the external current being constant [2,15]; it will be 
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expanded in powers of h rather than the coupling constant g and depends 
on the new parameter X a , where a = 1, 2, 3, namely the path average, 



V eS (X a ) = Y^ nlV{l) (Xa) 



1=0 

= V{X a ) + ^Y,^G x \ + V^\x a ). (56) 

a 

where the superscript / indicates the number of loops involved in Feynman 
diagrams, the trace-logarithm functions are given by the ground-state en- 
ergy of harmonic oscillators, they are connected to the partial frequencies 
by 

^trlnG- = % 

2 x a 2 



= ^d x V D (X a ). (57) 



The frequency of the propagator is given as the sum of all partial fre- 
quencies 



a 



^/( int ) (X a ) is called interaction potential [2,13-15] and contains all one- 
particle irreducible vacuum diagrams. 

5.1 Weak-coupling limit 

The aim of this section is to deduce weak-coupling ground-state energy for 
the potentials (2) and (3) by using the effective potential [10,12]. The com- 
putation will be performed until the third-loop order at low temperatures, 
i.e. (T->0). 

The frequency of the propagator of the potentials (2) and (3) are now 
given, using (57) and (58), by 

u 2D = ^2u 2 + 2igX. (59) 
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^3D = V&O- 

The corresponding interaction potentials are read 



V^ t] (X,Y) = -\imk B T 



1 

36 



+ 



1296 



+ 



V™ (X, Y, Z) = -Vmk B T 



1 1 

72 f 1728 



648 



+ 



5184 



(60) 



+ 0(h 4 ). (61) 



O (h*) . 
(62) 



where the previous diagrams are deduced from new-Feynman laws 
• for propagator 

-»■ G u (r a , r b ) . (63) 



for vertices 



—D 

IT 



D D D 



(64) 



a=i a=j a=k 



Substituting (59)- (62) into (56) and taking into account (63) and (64), 
we obtain the expressions of effective potential at low temperatures 

limK ( ff 2D) (X, y) = y (X 2 + y 2 ) + igXY* + \u> + ly/u* + 2igX+ 



2^2 



3 (2uo 2 + 2igXf /2 



2 1 22 5 
+ 



3 648 27 1296 



324 v / 2ftV 
(2uo 2 + 2igX) 9/2 



+ 0(h 4 ). (65) 



hmy e ( ff 3D) (X,Y,Z) 



T^ X +Y +Z ^ +igXYZ+ — + 2^ 



2 1 



22 5 



3 5184 271728 



3^,4 



WtTg 



O (h A ) . (66) 
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The path average X a in the case of complex potentials can be expanded 
in the form [13,15] 

X a = i (X a0 + hX al + h 2 X a2 + ---). (67) 

Inserting this identity into (65) and (66), and performing the first 
derivative with respect to X a and setting the resulting expressions to zero, 
we get 



• 2-D 



9 v _ v 5g 



3 



X\o = X = 0; Xu = Xi — — X 12 = X 2 — ^^( 68. a) 

X 20 = Y = 0; X 21 = Y 1 ^0; X 22 = Y 2 = 0. (68.b) 

• 3-D 

X ak = 0. (69) 

for a = 1, 2, 3, and k — 0, 1, 2. 

Re-inserting (68) and (69), respectively, into (65) and (66), and ex- 
panding until the third-order in K yields the ground-state energy at low 
temperatures 

uo) . 5 g 2 h 2 223 g 4 h 3 „ , fi . 



24 M 3 w 4 864 M 6 co 9 

r 0] 3 , 1 g 2 h 2 7 g 4 h 3 „ , BN 

^ = + dk?" 576^ (71) 
The results are in agreement with those obtained in sections 2 and 3. 



5.2 Strong-coupling limit 

Let us now resume the effective potential. Since the effective potential is 
expanded in powers of h, Kleinert's square-root will be converted accord- 
ingly to [13-15] 

uj -> ttVl + hr. (72) 

with 

r = "Mp- (73) 
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5.2.1 Two-dimensional case 

First-order Substituting (72) and (73) into (65) at first order in h, we 
obtain 

Kg? y) = y( x2 + y2 ) + % v xy2 + t + 5 v^T^x. (74) 

We now optimize in f2 and X a ; the resulting equations allow us to 
determine the strong-coupling behavior of X a given by [13,15] 

X a = -ig- lth (X a0 + X al g' 4 / 5 + X a2 g~^ + •••). (75) 

where the corresponding coefficients read as 

X = 0.417288 ff /5 ; X 1 = -0.245404 u 2 ; X 2 = 0.063624 <j 2 7r^6.a) 
Y k = V2X k ; fc = 0,l,2. (76.b) 

Re-inserting the results (76) into (75), and again in (74), we obtain the 
strong-coupling behavior of the ground-state energy 

4 1} ss 1.1263168 ff /5 . (77) 
Second order In second order, the effective potential becomes 



V%?(X,Y) = ^{X 2 + Y 2 ) + rgXY 2 + \ 



lo 2 „ u 2 + tt 2 + AigX 

— + + 

^ 2 ^Jtt 2 + 2igX 



a 2 h 2 

9_± (7g) 

12 {Q 2 + igX) 2 y ' 

Following the same steps as was done in the first order, the strong- 
coupling behavior of the ground-state energy in second order is 

« 1.13595605 ff /5 . (79) 
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5.2.2 Three-dimensional case 

First order Substituting (72) and (73) into (66) at first order in h, we 
obtain the effective potential 

Kg? (*) = y (X 2 + Y 2 + Z 2 ) + igXYZ + ^. (80) 
Let us now optimize in Q and X Q , with a = 1,2, 3. We can note that, 
from (80), the first derivative with respect to f2 is constant I = — I and 



2 f 

different than zero. In order to avoid this ambiguity, we need to substitute 
Q in (80) by a mathematical trick according to 



ft -> v 7 ^ 2 + 2 ^AX, (81) 

where A < 1 and by performing A-expansion around zero (A — > 0), we 
obtain 



AV 648' 1 5A' ' 100A 

y = 0; Y x = u 2 ; Y 2 = h' 6/5 . (82.b) 



1 5/ ^ 2 ^ 2 3cu 4 J 81 , . 

Z " = "V 648 ; Zl = "5A ; Z2 = "50aV5F- ^ 

Re-inserting (81) and (82) into (80), taking into account (75), we obtain 
the strong-coupling behavior of the ground-state energy at first-order 

Q3/5 

e W = __ ^6/5 w x 5702317 £6/5 _ ^ 83 ) 

Second order Substituting (72) and (73) into (66), and expanding the 
result until the second-order in h yields the effective potential 

3 (uj 2 + ft 2 ) h h 2 



V%? (X, Y,Z) = ^ (X 2 + Y 2 + Z 2 ) + igXYZ + 3( " 2 +/ )n + 



4ft 24ft 4 
(84) 
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Again, the strong-coupling behavior of the ground-state energy in sec- 
ond order is obtainable following the same steps as was done in the first 
order, we obtain 

4 2> = l^ /'V^ 5 fi6/5 * h5783441 ft6/5 - < 85 > 

In both cases, one can sees that the new values of the leading strong- 
coupling coefficients are in good agreement with those expected with a 



small deviation from the first order, i.e. 
0.5%, respectively. 



£ (2) _ (1) 
fc fc 



,(1) 



, in order of 0.8% and 



6 Conclusion 

In sections 2 and 3, we have derived the weak-coupling coefficients for 
the ground-state energy for the potentials (2) and (3) using the Feynman 
diagrammatical expansion and Bender- Wu recursion relations. Both re- 
sults are in agreement with those obtained by the old-fashion perturbation 
theory, namely the Rayleigh-Schrodinger method. We have proceeded in 
section 4 to resume the weak-coupling series (70) and (71) to the strong- 
coupling limit by applying Kleinert's square-root trick. However, the lead- 
ing strong-coupling coefficients for the ground-state energy lie below those 
calculated in weak-coupling limit, thus, the rate of convergence is less sat- 
isfactory. We introduce in section 5 the effective potential which, once 
combined with the VPT, allows to determine the same weak-coupling se- 
ries obtained in sections 2 and 3. In order to recover the rate of convergence 
in the strong-coupling limit, we proceed to resume the effective potential. 
It turns out that VPT-results for leading strong-coupling coefficients ob- 
tained in second order (79) and (85) approaches those expected with a very 
small deviation compared to the first order (77) and (83), thanks to the 
introduction of the path average variational parameter. 

Following [16], it is interesting to apply the covariant effective potential 
for particle moving in one-dimensional complex cubic potential for solving 
the corresponding Schrodinger equation with position-dependant mass. 
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